Cluster X-ray Gas 


Galaxies only contain about 20% of the baryons present in rich 
clusters. The other ~ 80% resides in the intracluster medium. 
These baryons can find their way into the intracluster medium in 
3 Ways: 


They can fall into the cluster as it forms 

They can be blown out of small galaxies by winds 

They can be swept out of larger galaxies by ram pressure 
They can be ejected from galaxies by the jets and winds of AGN 


Let’s consider the first three of these mechanisms one at a time. 
The first is simple. A rich cluster contains M ~ 10!°Mo of ma- 
terial within a ~ 1 Mpc radius. If an atom of mass my falls into 
this cluster from infinity, then its thermalization temperature is 
given by 


GMmyH 
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In the absense of cooling, this gas will emit x-rays. 


Galactic Winds 
[Mathews & Baker 1971, Ap.J., 170, 241] 


Consider an elliptical (or spheroidal) galaxy of luminosity, £, and 
mass-to-light ratio, Y. The stars in the galaxy are losing mass at 
the normal rate of Mrsm ~ 10-1! Mo per year per unit solar lu- 
minosity of the galaxy. Because the stars are moving isotropically 
within this galaxy, the gas will thermalize at a temperature given 
by 


GYLm 1 3 
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In other words, 
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Now assume that supernovae occasionally go off in the galaxy, and 
eject Magn per year per unit solar luminosity into the ISM/null. 
This matter will likely be thrown off at ugy ~ 10,000 km s7! 
and thermalize with an efficiency, e. (In other words, 1 — e of 
the bulk motion energy of the ejecta will be radiated away.) The 
temperature of this ejecta, if thermalized, would be 
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When this (much hotter) material mixes with the ambient ISM, it 
heats it. The final temperature of the system will become 


i T ANT. 
ER Mism ISM + Msn SN (10.06) 
Misa + MsN 


2 


Explicitly, this is 
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Next, we consider the escape velocity for the system, 
2GM 
esc — tS 10.08 
v z (10.08) 
which means that the escape temperature is given by 
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Now we ask the question: when is the ISM temperature greater 
than the escape temperature? 
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which, with a little math, becomes 
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For reasonable numbers, such as e€ ~ 0.1, Mrsm ~ 107 Mo, 
Men ~ 10715 Mo (both per year per solar luminosity), Y ~ 10, 
and R ~ 10 kpc, Lying ~ 2 x 101° Lo. Small galaxies therefore 
cannot hold onto their interstellar medium. Once a supernova goes 
off, it will heat the ISM, and create a galactic wind (sometimes 
called a Mathews-Baker wind). 





Galaxy Sweeping 
[Gunn & Gott 1972, Ap.J., 176, 1] 


Interstellar matter may be swept out of a galaxy by ram pressure, 
as the classic calculation of Gunn & Gott shows. To repeat their 
analysis, consider a spiral galaxy moving through a cluster that 
already has a low-density x-ray gas. Also, following their prescrip- 
tion, assume that most of the mass of the spiral resides in its disk. 
(This is not true, of course, but since the paper was written before 
the discovery of dark matter halos, we’ll go along with it.) 


Now let’s approximate a spi- 
ral galaxy as an infinite plane, 
with surface mass density o,. h 
From freshman physics, the 
force required to lift a particle 
of mass m from the surface of 
the galaxy to infinity is 


ie 
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Now let the surface density of gas in the plane of the spiral galaxy 
be Ogas. Since the mass of gas contained within a surface area 
dA is related to the surface density by dm = Ogasd A, the pressure 
required to strip gas from the galaxy is 





dm 
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Now when a galaxy moves through a cluster, the gas in the disk 
will feel a ram pressure due to x-ray gas. If the density of the x-ray 
gas is pz, then, from fluid mechanics, the ram pressure is 


P = prv’ (10.14) 





Consequently, the gas inside the galaxy will get stripped out by 
the ram pressure when 





Prt’ > NGOO gas (10.15) 


This is the Gunn-Gott criterion for sweeping of gas from spirals. 
When does this occur? Well, there are 2 x 10"! stars in the Milky 
Way, and the size of the Milky Ways disk is R ~ 20 kpc. Ifa typical 
Milky Way star has one solar mass, then the approximate stellar 
surface density is 0, ~ 2 x 10!!/7R? ~ 0.03 gm cm™?. Also, if the 
density of particles in the local interstellar medium is ~ 1 atom 
cm? and the thickness of the disk is ~ 200 pc, then the surface 
density of gas is Ggas ~ 6 x 10% atoms cm~*. Plugging these 
numbers into (10.15), and picking a typical velocity for a galaxy 
moving through a cluster (v ~ 1300 km s~'), yields a requisite 


x-ray gas density of ~ 5 x 1074 atoms cm™?. 


This calculation underestimates (by up to an order of magnitude), 
the density needed strip a galaxy, since it ignores the presence of 
dark matter. However, the qualitative result is the correct: it does 
not take much of an intercluster medium to strip spiral galaxies 
of their gas. When this occurs, star formation in the spirals will 
rapidly come to a halt (since they’ll be no gas to form stars), and 
the galaxy will (presumably) change into a lenticular. At the same 
time, the increased amount of intracluster material will cause the 
next spiral galaxy to be stripped that much easier. As time goes 
on, the amount (and x-ray luminosity) of the intracluster medium 
can be increased substantially. 





The X-ray Gas 


What will be the state of the x-ray gas? Given that the galaxy 
velocity dispersion in a cluster is ~ 2000 km s~+, a typical virial 
temperature of the gas will be T ~ 108 K. Also, for reference, a 
typical density for intracluster gas is ne ~ 107? cm7?. 


TIMESCALE FOR EQUIPARTITION 


Let us first ask how long it takes an electron (or proton) to share 
its energy with its surrounding and come to thermal equilbrium. 
It is fairly straightforward to show that the mean free path of a 
charged particle in a plasma is 


33/2 (kT)? 
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where e is the charge of the electron (or proton), n is the number 
density of particles and A is related to the Debeye length, i.e., the 
ratio of the largest to smallest impact parameters. (The derivation 
of this equation is straightforward; you will undoubtably see it, or 
the equivalent for a gravitational plasma of stars, in another class). 
Numerically, A is 


T Ne —1/2 


Note that this is a log quantity, hence it is virtually independent 
of ne and T. When we plug in the numbers, the mean free path 
of an electron (or ion) is 





A k T k k 10.1 
EAEE (=) T o ane) 


This is almost the length of a galaxy! 


Now, a reasonable estimate of the equipartition time is 


teq = /(Urms) = A/ ~ (10.19) 


(You can do the computation in a lot more detail, but the result 
will stay the same, to a factor of a couple.) Thus, from (10.16) 


= 3m1/2 (kT)? 


“1 —— Art/2netln A (1020) 


For proton-proton interactions (m = mp), this works out to be 


T aie n -1 
3 T Pp p 
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while for electron-electron interactions (m = me) 
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teq = 3.3 x 10 yr (5) (m) (10.22) 








Note that these timescales are very short (cosmologically speak- 
ing). Consequently, any new matter introduced into the intra- 
cluster medium will shortly come into equipartition with that 
medium. 


TIMESCALE FOR COOLING 


The next question to ask concerns the cooling timescale for the 
x-ray gas. In other words, how long with it take the temperature 
of the gas to decrease by 1/e. 


dT dinT\~? 
coo iw ee SS 10.2 
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Because most of the ions in an x-ray gas are either hydrogen- 
or helium-like, there are no low-lying levels to enable collisional 
cooling to work. Thus, the principal cooling mechanism is thermal 
Bremstrahlung. Recall that energy emitted by free-free emission 
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Thus, 
Eff = = = snk = 10a, T)/? ergs cm~’ s7! (10.25) 
which means that 
1/2 1/2 —1 
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This is more than a Hubble time. So, unless we’re talking about 
gas in the very center of the cluster (where the density may be 
high), the x-ray gas will not be cooling significantly. 





TIMESCALE FOR PRESSURE WAVE PROPAGATION 


Finally, let’s ask how long it takes a pressure wave to cross the 
cluster (and distribution any inhomogeneity over the entire sys- 
tem). This number is simply the size of the cluster divided by the 


sound speed, t.e., 
| ykT 
taS by ore =D — 10.27 
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where y is the ratio of specific heats (y = 5/3) and uy is the mean 
molecular weight (of the order of one). Plugging in the numbers 


gives 
GN ae T) 
= 8 2 cee 10.2 
tp = 6.5 x 10 (=) (=) (10.28) 


Once again, this is short compared to a Hubble time. Thus, the 
gas will be in pressure equilibrium with its surroundings. In other 
words, hydrostatic equilibrium will hold. 





X-ray Gas and Hydrostatic Equilibrium 
[Fabricant et al. 1980, Ap.J., 241, 552] 


Given that an x-ray gas is in thermal and pressure equilibrium, 
it is fairly straightforward to use the distribution of gas to derive 
the mass distribution of the cluster. Begin with the equation of 
hydrostatic equilibrium 


dP | SET 





aie ee 10.2 
dr r2 (10:23) 
and the ideal gas equation 
Poc pr (10.30) 


HMH 
If you take the derivative of (10.30) with respect to r, then 


dP k dp dT\  GM(r) 
= wt $ 0} =— p (10.31) 





2 
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Now if we multiply through by r/p, and take T outside the bracket, 
then 





dP kT {io rE, =- 


dr umy pdr T dr r 


= kT(r) [dnp dT) Mr) 
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= kT(r) [dnp dlnT 
MAS — umygG ie s dlnr b ee) 














So 











Note that if the gas is truly isothermal (due to the short equipar- 
tition timescale), the last term in the bracket disappears, making 
the mass a function of only one derivative. 
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Converting from Projected Radius to True Radius 


Unfortunately, equation (10.33) is of limited use. In particular, the 
density, p, is not an observable quantity; what is measured is the 
projected x-ray emissivity, [,(R), on the plane of the sky. How- 
ever, it is possible go from projected intensity to three-dimensional 
emissivity, €, (r) via an Abel integral. 


If the emitting x-ray plasma is spherically symmetric, then the 
flux one observes at any projected position in the cluster is 


l > 





L,(R) =2 in e (r)dz = 2 L a (10.34) 


R Be eS 


Now let’s introduce a dummy variable € (which takes the place of 
r), and multiple each side of the equation so that 


RdR œ  e,(r)rdr RdR 
I (R) (R2 — ¢2)1/2 = i (r2 — R2)1/2 (R2 — €2)1/2 (10.35) 
Now let’s integrate both sides over dR from € to infinity 
~° I, (R)RdR ee 
| (R2 — €2) (R2 _ £2\1/2 of rar | = a ate — €2)1/2 
(10.36) 


Now comes the tricky part. Let’s switch the order of the inte- 
gration. As the figure below shows, in order to do this, we must 
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also change the limits of the integration. In the first case, r is 
integrated from R to infinity (parallel to the x-axis) and R is in- 
tegrated from € to infinity. In the second case, to cover the same 
area, R goes from € to r, and R is taken from € to infinity. Thus 


°° 1,(R)RdR _9 Gi d . RdR 
¿ (R? — £2) 1/2 — é P e (r? — R?)1/2( R2 — £2) 1/2 


(10.37) 
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Now the second integral on the right hand side can be evaluated. 


If we substitute ; ; 
De ae 


e (10.38) 


For the integral is over R, r is held constant, so 


RdR = (r? — €?)ndn 
| (r2 — R2)1/2(R2 — €2)1/2 | (r2 — €2)1/2(1 — 2) 1/2(r?2 — €2)1/2 


1 
dn T 

= — = 10.39 
|, Te? oa 

ee © T (R)RdR oo 
‘ TR? EA = | e,(r)rdr (10.40) 

or, if we differentiate, 
ld [{* RdR 

sa a E 10.41 
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Converting from Density to Flux 


The projected x-ray flux distribution can be converted to the 
true x-ray flux distribution via the Abel integral. To go from 
flux to density, consider that the x-ray emissivity from thermal 
bremsstrahlung is 








32re? ( In \1? z 
Z z2 eNi —hv/kT4q = 2A T 
Eff 3mec’ (r) Nen i Gf fe V Me ( ) 
(10.42) 


Consequently, 


(ere) = (n,A(T)) = p°(A(T)) (10.43) 





where the angle brackets represent an integral over a energy band- 
pass (from vı to v2). Thus 





din(erp) _ Iln p dla (A(T) 


dlar dlnr dlnr ae) 


Now we solve for p and re-write the derivative of A(T), so 


dln p _ ld(err) _ 1din(A(T)) dln T 


dlor 2 dlnr 2 dinT dilar (Oe) 











If we then substitute this in (10.33) we get 











TES kTr fdlaT 1dln(eşp) ~ A din(A(T)) dlnT 
~amaG \dinr 2 dinr 2 dint | doir 
(10.46) 


or 





_pmAGM(r) S 1dln{esp) = dla T Pe 1dln(A(T)) 
kTr 2 dlnr dlnr 2 dlaT 
(10.47) 
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Note the terms in the above equation. The derivative of A(T) 
details how the free-free emission in an x-ray band changes with 
temperature; it is well known from basic physics via (10.42). Also 
known is the derivative of the x-ray emission with radius; this 
comes from the Abel inversion of the observed surface brightness. 
Therefore, given a function for T(r) and a boundary condition 
(usually the central temperature or pressure at infinity), this or- 
dinary differential equation can be solved for mass. 
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Isothermal X-ray Gas 


If an x-ray gas is in hydrostatic equilibrium, then 


dP GM(r) _ dQ 


dr r2 as. 





(10.48) 


where Q is the local potential. In addition, if we take the derivative 
of the ideal gas law, 





dP kT dp 
— = —— — 10.49 
dr umy dr ( ) 
Putting these two equations together yield 
kT dl dQ) 
cae cecal lee (10.50) 





umy dr dr 


Now consider: the x-ray gas isn’t the only thing effected by the 
cluster potential. The galaxies are also moving about. Let’s define 
a characteristic “temperature” of the galaxies, Tal, via 


1 3 
5 bmn l’) = 5kTeal (10.51) 


where (v?) is the mean square velocity of the galaxies. Note, how- 
ever, that this is not the observed velocity dispersion: since we 
only measure radial velocities, all we can see is one component of 
the motion. The observed velocity dispersion, @gai, is related to 
the true dispersion by (v?) = 3¢ga1. (Of course, this is only true 
if the orbits are isotropic ((vz) = (vy) = (vz)) and the cluster is 
virialized.) So, from (10.51) 


2 
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1 3 
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If the galaxies are in virial equilibrium, then this temperature is 
related to the potential by 


uma 

3k 
But how does this temperature compare to that of the x-ray gas. 
One can imagine that the x-ray gas may be cooler than the galax- 
ies, due to its Bremstralung radiation, or the inclusion of non- 
virialized, possibly foreground or background galaxies in the ve- 
locity dispersion measurement. Or, one can hypothesize that the 
gas is hotter, due to infall. (Recall that the kinematic energy an 
object has from the virial theorum is half that of a particle dropped 
from infinity.) 


2K+N=0 = —= Tal = Q (10.53) 


2 
Osal Z F 
ga 3 


To parameterize the problem, let’s define @ as the ratio of the 
virial temperature of the galaxies to the x-ray temperature of the 
gas, t.e., 





LNA 
B= op Fal 
Since both the gas and the galaxies must feel the same potential, 


then, from (10.50) and (10.52) 





(10.54) 


T dl In Pga Q 
kT dln pz > dMpga_ d (10.55) 








umy d gal dr dr 


where pz is the density of the x-ray gas, and pga) the density of 
the galaxies. Re-writing this, we get 








dinp, uMH > dln pea d1n Peal 
BEP EA sG 10.56 
dr m ge O dr CORN 
If we trivially integrate over r, then 
Px = Phol (10.57) 
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The Total Mass of an X-ray Cluster 


If the gas is isothermal, then its radial distribution should be that 
of an isothermal sphere. In that case, a reasonable approximation 
for the density distribution over the inner ~ 2 core radii is 


( )= PO 
a= (1 + #2)3/2 





(10.58) 


where po is the central density, 7 is r/ro, and ro is the core radius. 
This is sometimes called a modified Hubble Law. 


The modified Hubble Law departs from an isothermal sphere at 
large radii. Nevertheless, it is common to extrapolate this law to 
infinite radius, and derive the total x-ray gas mass of the cluster 
in terms of the core radius and central density. (But note: this 
only works if @ > 1; otherwise, the mass is infinite, even with this 
approximation.) The calculation is fairly straightforward. One can 
measure ro by counting galaxies and fitting the radial distribution 
to an isothermal profile (or a Hubble law). One can also measure 
the central density from the observed amount of x-ray emission. 


So 
Mz = [am z eee, 
(1+ (r/r9)2)2°? 


= anpor f P2 (147P) ap (10.59) 
(0) 
If we let z = 77, then 


Fa = 
Mz =4npor§ | 50 (1+ zx) 3/2 dg 
0 





= 2nporg 


3(G—1) 
se 2 ) 


2 AY 
Poro r (2) 


= 7°! (10.60) 


or 


w f n VTE) 
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where T is the gamma function. 
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6 and the Flux Distribution of an X-ray Cluster 


One can determine O by directly measuring the temperature of 
the intracluster medium and comparing it to the velocity disper- 
sion measured for the galaxies. However, an alternative method, 
which is often used is to integrate (10.34) over the line-of-sight, 
and compare the result to the observed x-ray surface brightness. 
From (10.34) this surface brightness should be 


oo 2 
A P 
rR (r2— R) 
(10.62) 
If we substitute for n? using (10.57) and (10.58), and let z = 


(r? — R?)/(r2 + R?), then 
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When one fits the x-ray surface brightness for many clusters in 
this manner, one finds (3) ~ 0.65. On the other hand, when 
one measures the galactic velocity dispersion and compares that 
number to the gas temperature as determined directly by fitting 
x-ray spectra, then (3) ~ 1.2. Either 


a) The galaxies and/or the hot gas are not isothermal. 
b) The galaxies are not in isotropic orbits. 


c) The galactic velocity dispersion is overestimated by the inclu- 
sion of contaminating foreground and background galaxies. 


d) The x-ray gas and galaxies do not see the same potential. 
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Fig. 21. Upper panel: mean temperature profiles obtained from 

this work (black circles, LM08), by De Grandi & Molendi (blue 

0.6 | s | : | : squares, DMO02), by Vikhlinin et al. (red upward triangles, V05), 

4 6 8 10 and by Pratt et al. (green diamonds, P07). All profiles are 

Temperature (kev) rescaled by k7m and Rjgo as defined in Sect. 4, The dashed 

line shows the best fit with a linear model beyond 0.2 Rigo (see 

Fig. 10. 8-temperature relation. £ is the parameter obtained when fit- Sect. 6.1) and is drawn to guide the eye. Lower panel: residuals 

ting with a 8—model the X-ray surface brightness profile. There is a with respect to the linear model. The LM08 profile is the flattest 
trend to find larger £ values for larger temperatures. one. 


[Castillo-Morales & Schindler 2003] [Leccardi & Molendi 2008] 











20 


